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We provide for the first time accurate assessments of the consequences of violations of self- 
consistency in the Hartree-Fock based random phase approximation (RPA) as commonly used to 



causes a spurious enhancement of the breathing mode energy for spin unsaturated systems. Con- 



^ ' calculate the energy Ec of the nuclear breathing mode. Using several Skyrme interactions we find 

cn ■ 

that the self-consistency violated by ignoring the spin-orbit interaction in the RPA calculation 

o 

^—N I trarily, neglecting the Coulomb interaction in the RPA or performing the RPA calculations in the 

TJ scheme underestimates the breathing mode energy. Surprisingly, our results for the ^^Zr and 

r^ , ^°*Pb nuclei for several Skyrme type effective nucleon-nucleon interactions having a wide range 

^ , of nuclear matter incompressibility {Knm ~ 215 — 275 MeV) and symmetry energy ( J ~ 27 — 37 

MeV) indicate that the net uncertainty {5Ec ~ 0.3 MeV) is comparable to the experimental one. 

, 5^ ; PACS numbers: 21.60.Jz, 24.30. Cz, 21.65.+f 



The Hartree-Fock (HF) based random phase approximation (RPA) provides a microscopic 
description of the nuclear compressional modes. The most special of these compressional 
modes is the isoscalar giant monopole resonance (ISGMR) also referred to as the breathing 
mode. The centroid energy E^ of the ISGMR enables one to determine the value of nuclear 
matter incompressibility coefficient Knm which plays an important role in understanding a 
wide variety of phenomena ranging from heavy-ion collision to supernova explosions. Recent 
experimental data \v\ for the Ec in heavy nuclei have rather small uncertainties (~ 0.1 — 
0.3 MeV). Since, the uncertainty 6Ec associated with E^. is approximately related to the 
uncertainty 6 Knm in Knm by, 

K ~ E ' ^ ' 

the value of 5Knm is only about 10 MeV, for Knm = 250 MeV and E^ = 14.17 ± 0.28 

MeV for the ^''^Pb nucleus. On the theoretical side, most of the results obtained within 
the HF-RPA theory, as employed for the determination of Ec, are plagued by the lack of 
self-consistency, see however, Ref. 2|- Self-consistency is violated due to the neglect of the 
spin-orbit and Coulomb terms in the particle-hole interaction used in the RPA calculations. 
Furthermore, some of the RPA calculations are performed in the TJ (isospin) scheme. Un- 
certainties in the value of Ec calculated in the HF-RPA approach arise also due to various 
numerical approximations such as the use of very small particle-hole excitation energies and 
;he introduction of the smearing parameter F to smoothen the strength function. In Ref. 
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we have investigated in detail the effects of these numerical approximations and empha- 
sized the need to perform a highly accurate calculation of the strength function. So, it is 
of utmost importance that in order to determine an accurate value of Knm one must have 
an accurate knowledge of the effects of violations of self-consistency on the centroid energy 
of the ISGMR. We may also point out that modifying the particle-hole interaction in an ad 
hoc manner in such a way that the spurious state associated with the center of mass motion 
appears at zero energy |^ does not restore the self-consistency (see also Ref. [3|). 

In this work we provide for the first time accurate assessments of the effects of violations 
of self-consistency on the constrained and the scaling energies. 



Econ= \l — -, and E,= .—, (2) 
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of the ISGMR. Here 
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rrik = 00 S{uo)duj (3) 

Jo 



is the energy moments of the strength function 

Siuj) = J2\{^\f\^)\'Si^-^n), (4) 

n 

for the monopole operator /(r) = Yli=i ^I- 

In the following the fully self-consistent and highly accurate values of the constraint 
energy, Econ, and the scaling energy, Es, are obtained by calculating the energy moments 
ruk for k = —1,1 and 3 using the constrained HF (CHF), the double commutator and the 



generalized scaling approaches, respectively 



5| . The self-consistent values of Econ and Eg are 

in order to assess the 



then compared with those obtained within the HF-RPA approach 
effects of each type of violation of self-consistency. It may be pointed out that very recently 
in Ref. j6| some preliminary attempts are made to estimate the effects of self-consistency 
violation on the value of Econ- However, Eg was not considered and the different components 
contributing to the self-consistency violations are not dealt separately. We show here that 
the appropriate knowledge of the effects of each of these components on the ISGMR energies 
are instrumental in validating the HF-RPA calculations. 
We consider here the HF solution for the Hamiltonian, 

^"Er. + jEf'i + jE'^f. (5) 

where T is the kinetic energy operator, V'" is the Coulomb interaction and V is the effective 
two-body interaction of the Skyrme type [3], 



+ ^ti(l + XiPf2)x 



k ?2(5(ri - ra) + (5(ri - rs) A; ^ 
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+t2 (1 + X2P^2) k i2'5(ri - rs) k 12 

+ ^^3 (1 + XsPf^) P" (^4^) '^(ri - r2) 

+iWoki26{ri-r2){o\+a^) X ku, (6) 

where P{^ is the spin exchange operator, a « is the Pauli spin operator. A; 12 = — i( Vi — V2)/2 
and k 12 = — ^(Vi — V2)/2. Here, the right and left arrows indicate that the momentum 
operators act on the right and on the left, respectively. The parameters of the Skyrme force 
ti, Xi, a and Wq are obtained by fitting the HF results to the experimental data for the bulk 
properties of finite nuclei. 



Once the HF equations are solved, the mi and rris can be expressed in terms of the 
ground state density p, kinetic energy density r and spin density J. The value of mi can be 



obtained from the corresponding double commutator 



and is given as, 



mi = 2 — (r ), 
m 



(7) 



where, 

(r^) = / r^p{r)dr. (8) 

The cubic moment m,-^ can be evaluated in the generalized scaling approach and it is given 
by 0,3, 

1 f2h^V 
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(3a + 2)(3a + 3)Ep], 
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Efin = / [aopr + aipiTi + c{y pf + diy pif] dr 
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1,3,- _i_,3(,, + !),«,; 



dr 



Eso = -^Wo / [pV • J + PpV ■ Jp + p„V ■ J„] dr 
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with pi = pn- pp, n = Tn - Tp, ao = jgiSti + (5 + 4^2)^2), ai = 1^2 (a;2 + i) - ii(a;i + |), 
c = ^(9ti — (5 + 4x2) + ^2) and c? = — ^(3ti(a;i + |) + t2(a;2 + ^)). As described in detail in 
Ref. |5|, m_i can be evaluated via the CHF approach and is given as. 
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where, (Ox) = ($a \0\ $a), with $a being the HF solution for the constrained Hamiltonian, 
Hx = H- A/. 

The corresponding HF-RPA values of mk are obtained from Eq. (jH)) using. 



S{uo) = -Im[Tr{fG{uj)f)] 



(16) 



with G{uj) being the RPA Green's function [21, y|. The differences between the values of nik 
and consequently the difference 

6E = E{RPA)-E{SC), (17) 

between the energies obtained within the HF-RPA and the self-consistent (SC) values ob- 
tained within the HF approach are due to the violation of self-consistency in the HF-RPA. 
Since we are unable to obtain a fully self-consistent value for mo, we only provide the value 
of the centroid energy 

E. = ^. (18) 

mo 

obtained within the HF-RPA. Thus we use SEcon and 6Es as measures of 6Ec. We point out 
that we find Ec ~ Econ- 

Before we present the main results, some remarks about the numerical accuracies are in 
order. We carry out the HF calculations using a box of 18 fm and a mesh size of 6r = 0.02 
fm. To determine r7i_i given by Eq. ()15j) with an accuracy of 0.1% we calculate the first 
derivative of (rj) using a five-point formula with the increment 6X = 0.02. In order to 
limit the spurious enhancement in rrik ( particularly m^), due to the Lorentzian smearing 
of the RPA Green's function, to below 0.2%, we carry out the continuum RPA (CRPA) 
and the discretized RPA (DRPA) calculations using the very small smearing parameter of 
r/2 = 5 X 10"^ MeV and the maximum energy of AE^. in the integral of Eq. Q. In Ref. [3| 
we show that in order to calculate the centroid energy accurate to within 0.1 MeV, one must 
set the cut-off for the p-h excitations -E^""^ > 400 MeV. In what follows, we shall present 
our DRPA results obtained with E^"^ = 950 MeV. 

In Table |l] we demonstrate the accuracy of our calculations by employing a simplified 
Skyrme interaction ^ with to = -1800 MeVfm^, tg = 12871 MeVfm^ and a = 1/3. We 
consider the ^°Zr nucleus and keep the spin-orbit and the Coulomb interactions switched off 
in the HF as well as in RPA calculations, so that the HF-RPA calculations are fully self- 
consistent. Self-consistent (SC) values of rrik and the energies Econ and Eg shown in the first 
row of this table are obtained using Eqs. (0), ©, © and (fT3j). These values are accurate 
to within 0.1%. We compare these highly accurate results with the corresponding ones 
obtained within the HF-RPA approach. The CRPA and DRPA results presented in Table 
m are calculated by integrating the RPA strength function appearing in Eq. Q up to the 
different values of the maximum energy Umax ■ One can easily verify that the CRPA as well as 



DRPA results are quite close to the highly accurate SC results (first row). For instance, the 
maximum uncertainty is associated with the Eg and it is about 0.3% for uJmax = 100 MeV. 
Since, the RPA energies Econ, E^ and Eg obtained by integrating Eq. Q up to ujmax = 100 
and 120 MeV are very close, we shall adopt in the following the value of (jJmax ~ 4£^c for the 
RPA calculations. 

Having demonstrated in Table Q] the high accuracy of our numerical calculations, we now 
consider the effects of the violations of self-consistency. To investigate effects of each of the 
components contributing to the self-consistency violations, we carry out the calculations for 
some specific nuclei using several interactions [lOl, lUl llj, ll3| • We denote by 6is{6c) = or 
1 that the spin-orbit (Coulomb) interaction is turned off or on in the HF calculations. In 
the RPA calculations the p-h spin-orbit and Coulomb interactions are always absent. Thus, 
the HF-RPA calculations for symmetric nuclei (N = Z) with 6is = 6c = are fully self- 
consistent. The effects of ignoring the p-h spin-orbit interaction on the ISGMR energies can 
be singled out by carrying out the calculations for symmetric nuclei with 6is = 1 and 6c = 0. 
Similarly the effects due to the neglect of p-h Coulomb interaction can be estimated by 
performing the calculations for the cases with 6is = and 6c = 1. The HF-RPA calculations 
for asymmetric nuclei with 6is = 6c = yields the effect of working in the TJ scheme. In 
Table ITTl we present the results obtained using the SGII interaction [lO] for several nuclei for 
different combinations of 6is and 6c- Results for other Skyrme interactions exhibit similar 
features. The ^°'^°Ca and 80,1102^. nuclei considered here are spin saturated systems. But, 
^^Ni and ^°°Sn are spin unsaturated nuclei having I/7/2 and 1(79/2 orbits as the last occupied 
ones, respectively. As seen in Table HI it is once again evident that the fully self-consistent 
HF-RPA calculations for the ^°Ca and ^°Zr nuclei {6is = 6c = case) yield small values 
for 6Econ and 6Es (< 0.1 MeV). Considering the results for the cases with 6is = 1, we find 
that for the spin saturated nuclei the ISGMR energy do not get altered much even if the 
p-h spin-orbit interaction is not included. However, in case of the spin unsaturated nuclei, 
the effect of the neglect of the p-h spin-orbit interaction in the RPA leads to a significant 
spurious enhancement of the ISGMR energy. Following Eq. (^, one can easily verify that 
such enhancements in the ISGMR energies are equivalent to overestimating the value of 
Knm by about 30-40 MeV, which is unreasonably large in view of the current accuracy of 
the experimental data. Let us now focus on the results obtained for ^°Ca and ^^''Zr nuclei 
with 6is = 6c = and ^°Ca and ^°Zr nuclei with 6is = but 6c = 1. Compilation of these 



results enables us to investigate the effects of working in the T J scheme as well as the neglect 
of the Coulomb term in the p-h interaction used in the RPA calculations. One can conclude 
from Table |n] that both of these factors, contributing to the violation of self-consistency, 
tend to underestimate the ISGMR energy compared with its SC value. 

So far we have seen in a systematic manner how the inconsistency in implementing the 
HF-RPA approach introduces uncertainty in the ISGMR energy. Considering the results 
of Tables m one may expect that the total effects of the self-consistency violation on the 
ISGMR energy may be smaller. Because, the absence of p-h spin-orbit interaction tends to 
increase the ISGMR energy while working in TJ scheme or ignoring the Coulomb term in p-h 
interaction used in the RPA tends to underestimate the ISGMR energy. As an illustration, 
we present in Table ITTTl the results for ^^Zr and ^°^Pb nuclei which are widely used to extract 
the value of Knm from the ISGMR energies. These calculations are carried out for several 
Skyrme interactions such as SGII, SkM*, SLy4, SK255 and SK272 [liJ, Lll, UJ, M • The choice 



of these interactions cover a wide range of incompressibility coefficient Knm = 214—272 MeV, 
symmetry energy coefficient J = 27 — 37 MeV and the spin-orbit strength Wq = 95 — 130 
MeVfm^. We see that in the case of ^^Zt, the ISGMR energy Econ obtained in DRPA agrees 
with the corresponding self-consistent values to within 0.3 MeV. For the ^osp]^ nucleus, we 
have I 6Econ |< 0.4 MeV. For both nuclei, the behaviour of 6Es is similar to that of SEcon 
except for the fact that former one is larger by about 0.2 MeV. Thus, the net uncertainties 
in the ISGMR energies due to the self-consistency violations are reasonably small for the 
^°Zr and ^osp]^ nuclei and are comparable to that attained in recent experiments. Of course, 
this may not be the case for spin saturated nuclei near the drip line; the inconsistency in the 
HF-RPA calculations introduced by the Coulomb interaction and the asymmetry in these 
nuclei may overshadow the effect of neglecting the p-h spin-orbit interaction. 

In summary, we have investigated in detail the effects of the violations of the self- 
consistency in the HF-RPA calculation of the ISGMR energies Econ and Eg. In particular, 
we considered the self-consistency violations caused by ignoring the spin-orbit and Coulomb 
terms in the p-h interaction and by carrying out the RPA calculations in the TJ scheme. We 
performed the HF-RPA calculations for the ISGMR energies for several nuclei with the SGII 
Skyrme interaction and demonstrated that ignoring the spin-orbit term in the p-h interac- 
tion gives rise to a spurious enhancement in the values of Econ and Eg for spin unsaturated 
nuclei. On the contrary, neglect of the Coulomb term in the p-h interaction and performing 



the RPA calculations in the TJ scheme underestimates the ISGMR energies. Finally, we 
calculated the ISGMR energies for the ^°Zr and ^"^^Pb nuclei for the five different Skyrme 
interactions SGII, SkM*, SLy4, SK255 and SK272 and show that in these nuclei, widely 
used to extract the value of Knm.) the various elements contributing to the self-consistency 
violations tend to counterbalance their effects leading to an uncertainty of about 0.1 — 0.4 
MeV in the values of ISGMR energies, which is quite acceptable in view of the accuracy of 
the experimental data currently available. 

This work was supported in part by the US Department of Energy under grant 7^ DOE- 
FG03-93ER40773 and the National Science Foundation under grant #PHY-0355200. 
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TABLE I: Comparison of self-consistent (SC) values of the moments ruk (in MeV fm^) and the 
ISGMR energies (in MeV) for the ^''Zr nucleus calculated in the HF approach with the ones 
obtained within the continuum and discretized RPA using the smearing parameter T = 0.01 MeV. 
The RPA values are obtained by integrating Eq. ^ over the energy range — Umax (in MeV). 
The calculations are performed using a simplified Skyrme interaction with to = —1800 MeVfm^, 
ts = 12871 MeVfm"^ and a = 1/3. 





^max 


m_i 


mo 


mi 


mg 


^con 


Ec 


Es 


SC 




185.1 




96346 


54490227 


22.81 




23.78 


CRPA 


80 


185.0 


4189 


96229 


53549803 


22.80 


22.97 


23.59 


CRPA 


100 


185.1 


4190 


96289 


54024448 


22.81 


22.98 


23.69 


CRPA 


120 


185.1 


4190 


96303 


54182661 


22.81 


22.98 


23.72 


DRPA 


80 


185.4 


4194 


96359 


53721598 


22.80 


22.97 


23.61 


DRPA 


100 


185.4 


4196 


96426 


54292408 


22.80 


22.98 


23.73 


DRPA 


120 


185.4 


4196 


96449 


54544736 


22.81 


22.99 


23.78 



TABLE II: Comparison between the SC and DRPA results for the ISGMR energies (in MeV) 
for several nuclei obtained for the SGII interaction jlO|. The notation 9is{9c) = 1 or indicates 
that the HF calculations are performed with or without the inclusion of the spin-orbit (Coulomb) 
interaction. Note that the p-h interaction used in the RPA calculations does not include the 
spin-orbit and Coulomb terms. 



Nucleus 


9ls 


Oc 




DRPA 




SC 


5E,on 


5Es 


^con 


Ec 


Es 


^con 


Es 


40Ca 








22.53 


22.76 


23.85 


22.60 


23.92 


-0.07 


-0.07 




1 





22.31 


22.63 


23.96 


22.58 


23.93 


-0.27 


0.03 







1 


20.97 


21.24 


22.44 


21.38 


22.69 


-0.41 


-0.25 


802, 








19.57 


19.69 


20.26 


19.60 


20.32 


-0.03 


-0.06 




1 





19.49 


19.65 


20.36 


19.58 


20.33 


-0.09 


0.03 







1 


17.38 


17.54 


18.20 


17.81 


18.63 


-0.43 


-0.43 


56 Ni 


1 





23.14 


23.35 


24.44 


21.49 


22.41 


1.65 


2.03 


lOOSn 


1 





20.16 


20.27 


20.95 


19.04 


19.57 


1.12 


1.38 


60 Ca 








16.08 


16.81 


19.13 


17.32 


19.57 


-1.24 


-0.44 


ii02r 








15.64 


16.04 


17.35 


16.33 


17.52 


-0.69 


-0.17 



TABLE III: Comparison between the SC and the DRPA results for the ISGMR energies (in MeV) 
for the ^''Zr and ^"^^Pb nuclei. The DRPA results are obtained using the smearing parameter 
r = 0.01 MeV and the values of ujmax = 70 and 50 MeV for the ^*^Zr and ^osp^-, nuclei, respectively. 
Note that the values of SEcon Sive comparable to the experimental uncertainties. 



Nucleus 


Int. 




DRPA 




SC 


SEcon 


SEs 


i^con 


Ec 


Es 


i^con 


Es 


90 Zr 


SGII 


17.93 


18.05 


18.70 


17.70 


18.30 


0.23 


0.40 




SkM* 


17.72 


17.85 


18.55 


17.43 


18.00 


0.29 


0.55 




SLy4 


18.19 


18.33 


19.07 


17.86 


18.46 


0.33 


0.61 




SK255 


18.83 


18.97 


19.64 


18.87 


19.57 


- 0.04 


0.07 




SK272 


19.55 


19.70 


20.45 


19.49 


20.28 


0.06 


0.17 


208pb 


SGII 


13.53 


13.64 


14.22 


13.42 


13.94 


0.11 


0.28 




SkM* 


13.46 


13.61 


14.31 


13.36 


13.88 


0.10 


0.43 




SLy4 


14.04 


14.16 


14.80 


13.68 


14.22 


0.36 


0.58 




SK255 


13.83 


13.98 


14.65 


14.26 


14.88 


- 0.43 


-0.23 




SK272 


14.42 


14.58 


15.28 


14.73 


15.42 


- 0.31 


-0.14 



